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1. Introduction 


1.1. Ordination as a method of indirect gradient analysis 


The continuum model of community organization assumes that individual species popu- 
lation densities are distributed according to unimodal, gaussian functions, and by continuous 
replacement along a single environmental gradient, form a corresponding coenocline (Wnurr- 
TAKER 1978). Ordination techniques, such as Principal Component Analysis (PCA) and 
Reciprocal Averaging (RA), mathematically collapse a species by sample data matrix to 
a space of fewer dimensions. Indirect gradient analysis (sensu Gauen 1982) determines how 
well gradients of community composition, i.e. coenoclines, are represented by the abstract 
axes or components produced by ordination, and how well these axes or components are 
related to actual environmental gradients. While ordinations are useful in the context of 
delimiting relationships between communities of low beta diversity (sensu WHITTAKER 1978), 
and short portions of complex-environmental gradients, constraints of linearity and gradient 
“orthogonality” (sensu Gauen 1982) governing these eigenanalysis techniques may not be 
compatible with communities of high beta diversity. Swan (1970) recognized that coenoclines 
become increasingly convoluted in ordination space as the rate of species turnover, i.e. 
change in beta diversity, increases along an environmental gradient. The departure from 
linearity is spread among the second, third and higher level ordination axes as artifactual 
quadratic, cubic and higher order polynomial distortions of the first axis, respectively. 
KENDALL (1975) and Fasam (1977) termed the phenomenon the “arch” or “horseshoe” 
problem. 

Gradient “orthogonality” is the second problem facing community ecologists who wish 
to use ordination methodologies. However, PCA and RA may be rendered ineffective as 
trend-seeking methods for detecting relationships between species abundance data and 
underlying environmental gradients when the gradients are not independent of one another. 

The linear structure of the principal component model can be more closely adhered to if 
data partitioning is employed prior to extraction of eigenvalues and eigenvectors. Data 
partitioning may obviate use of oblique (partially correlated) principal component rotation, 
since the biological importance of non-orthogonal solutions is not always clear (Gauen 
1982). In an elegant application of the data-partitioning strategy to the fire history of sub- 
arctic woodland, Jounson (1981) separated temporal and spatial components in the vegeta- 
tion dynamics by fitting species abundance data to a predominant gradient of fire recurrence 
to produce a successional ordination. A complementary habitat ordination (uncorrelated 
with the first) was constructed from the residual variation not accounted for by regression 
along the temporal fire frequence gradient. 
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1.2. The-problem 


Ordination has been frequently applied to species-rich microarthropod communities 
primarily as a transformation to reduce the high dimensionality of the data. Although 
microarthropods are arranged primarily along a gradient of soil moisture (e.g. BoNNET 
et al. 1975; Topa & Tanno 1983), individual collembolan and acarine species are likely to 
be arranged along some other conspicuous macro-environmental gradient, such as vegeta- 
tional composition, that covarie with soil moisture status (e.g. BLACKITH & BLACKITH 1975; 
Curry 1977; Parr 1977; Appison 1980; Prat & Massoup 1980; Wautuy & LEBRUN 
1980; PTOWMAN 1981a, b; Poursin & Ponce 1982). Except within the simplest of com- 
munities (GODDARD 1979a, b; Boorn & Usuer 1984; WEST 1984), it has proved difficult. 
using ordination and indirect gradient analysis, to untangle the effects of soil moisture from 
more subtle influences of secondary gradients (e.g. soil pH, nutrient ion concentration) on 
soil faunal community structure (USER 1976). In addition, obvious arch distortions present 
in the reduced ordination space may have a profound effect on subsequent interpretation 
of a few concurrently, but coarsely sampled environmental variables; Acari and Collembola 
are typically distributed below the scale attained by bulk measurement in a complex mosaic 
of microhabitat patches. 

The purpose of the present study was to determine whether or not a data partitioning 
strategy could be successfully applied to gradient analysis and ordination of a species-poor 
collembolan community. The data for the analysis consisted of physico-chemical measure- 
ments and collembole species abundances which had been collected in the second field season 
(May to September 1980) of a 2-year study of soil faunal colonization of subalpine, reclaimed 
coal spoil. Data are also presented regarding the colonization of buried litter bag samples 
during the field season by microarthropods. 

The species abundance data were resolved into linear and non-linear components of varia- 
tion prior to ordination following the ideas presented by Jounson (1981). The first (non- 
linear) component relates species responses to a gradient of spoil moisture stress and its 
proxy variable spoil moisture content. This represented one aspect of the niches of the species 
involved in the ordination, i.e. extreme habitat adversity (sensu Souruwoop 1977). Residual 
variation not accounted for by moisture was ordinated, and represents species response to 
a gradient of potential spoil erosion. This represented a second aspect of the niches of species 
ordinated, denoted by Sournwoop (1977) as durational stability. 


2. Site description 


Three study plots (30 m «50 m) were established on bituminous coal overburden which had been 
seeded as pasture by the reclamation division of Cardinal River Coals Limited, Luscar, Alberta. The 
three sites differed in age since initial reclamation (8, 4, and 2 years). Five 2 m 10 m subplots were 
arranged, in turn, down the slope gradient within each plot. The ground cover was mostly vascular 
plant litter, underlain or supplanted by moss turf in poorly drained, low-lying areas of the plots. 
Surface runoff was intercepted by low spoil ridges (20 to 50 cm in height), spaced about 1 m apart, 
which ran discontinuously across the slopes, especially within the 2 year old plot. DuRALL ef al. 
(1985), and Parsons & Parkinson (1986) should be consulted for further desc riptions of the region, 
and reclamation history of the subalpine study sites. 

Each plot was surveyed in May 1979 to facilitate regular spacing of litter bags in an allied 
decomposition experiment (DenAUL etal. 1985; Parsons el al. 1986), and to provide topographic 
data for subsequent assessment of spoil erosional losses. The 1 mm mesh litter bags acted as crude 
sediment traps since burial of bags occurred with the onset of snow melt in 1980. Upslope angle, 
elevation (metres) above the toe of the slope, and aspect were recorded for each litter bag, and sub- 
sequently, for locations of spoil samples used in the present soil faunal study. 


3. Methods 
3.0. General 


Three core samples (9.8 em dia. « 5 cm) were collected at monthly intervals from each of the 
15 subplots, May to September 1980 inclusive. The procedures describing sample preparation and 
measurement of physicochemical variables, as well as extraction, enumeration and identification 
of Collembola are outlined in Parsons & Parkinson (1986). 
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3.1. Analysis of complex-environmental gradients 
3.1.1. Vegetation 


A simple analysis of mine spoil vegetational development was performed by ranking the moss 
cover and the litter cover in 0.25 m? quadrats located over the core samples collected in the present 
study. The ascending ranks (0—6) corresponded to increases both in percentage areal cover and in 
the thickness of the litter layer or moss mat over mineral spoil. Combinations of the ranks assigned 
to the litter cover and moss cover data formed broad, but distinct groupings; these “vegetation 
classes" were used in subsequent correlation and regression analysis of spoil water potential. The 
same vegetation classes were also used as “treatment” levels in non-parametric ANOVA of sediment 
density, pH, slope angle, and soil organic matter data, as well as categories in subsequent Multiple 
Discriminant Analysis of the ordinated collembolan data (Section 3.2.). Slope angles were converted 
to slope classes (Anonymous 1978) prior to ANOVA. 


3.1.2. Spoil water potential 


Spoil moisture stress was measured at monthly intervals using calibrated psychrometric thermo- 
couples (PT 51-10, Wescor Inc., Legan, Utah) attached to a dewpoint precision microvoltmeter 
(I R-33T, Wescor Inc.). The soil water potentials (y), expressed as positive values, were regressed 
against the ratio of actual volumetric spoil moisture content (Oy) to water content at saturation 
(Os). Vegetation class, pH, organic matter, and USDA textural classes of percent sand, percent silt, 
and percent clay were included in the regression analysis. Volumetrie moisture content at saturation 
was estimated as: 1 — (D»/Dp), where Dy is sample bulk density and Dp is mean spoil particle density 
(2.34 g cm-3). Spoil “water” or “moisture content" refers to the ratio [Oy/Os] in all subsequent 
analysis and discussion. All variables, except vegetation class and pH, were log, -transformed prior 
to regression analysis. Temporal differences in water potential were included in the regression ana- 
lysis by coding months as “dummy” variables (PEDHAzUR 1982). 


3.1.3. Spoil erosion 


Spoil erosional losses were predicted using equations obtained by regressing log, -transformed 
density of spoil (g cm~, O.D.W.) trapped in the 20 em x 20 em 1 mm mesh litter bags against tan- 
(slope angle), and log,9-transformed slope length (metres). Months in which the litter bags were sam- 
pled were coded as dummy variables. Mass losses of litter material in the bags and litter moisture 
content were rank correlated with actual sediment losses. Organic matter content, texture, moisture 
tension, and collembolan abundance and species composition of the spoil-infiltrated bags were also 
measured and correlated with sediment density. 

Predicted sediment density (SED) was rank correlated with spoil water potential, pH, organic 
matter content, and vegetation cover data derived from the spoil cores taken within the subplots. 
The pairwise correlations were tested for homogeneity among months by computing Chi-square 
statistics on the normalized coefficients (SoKAL & Ronurr 1969). 


3.2. Gradient analysis and ordination 


Several attempts were made to relate variation in collembolan density to variation in physico- 
chemical properties of Luscar spoil. In a direct gradient analysis (sensu WHITTAKER 1978), weighted 
averages of the untransformed collembolan abundance data were computed for each sample core, each 
month (May to September 1980). Ranks or weights averages assigned to each species corresponded 
to the order of increasing maximum species abundance along a gradient of increasing spoil moisture 
content as determined from data collected in the previous field season. During each month, samples 
were ordered as: qi = [Xi; * W;]/ Xi; where qu is the sample score for the ith core, Wj is the arbitrary 
rank given the jth species, Xij is the abundance of the jth species in the ith core, and XI. is the total 
collembolan abundance of the ith core. The **q-scores" were rank correlated with the environmental 
variables. 

Collembolan abundance data from the litter, 0 to 2.5 cm, and 2.5 to 5.0 cm strata were pooled, 
and transformed (logio X + 1). Species occurring with a frequency of < 4% within each set of 
45 monthly samples were excluded from trend surface and ordination analysis. The data remaining 
after the screen were standardized (by species) to normal scores. 

Trend surface analysis (TSA) extends the least-squares regression criterion for fitting a line to 
data in two dimensions to the fitting or mapping of a surface in three dimensions (Davis 1973; Ux- 
WII 1975). The simplest mapping function, a planar response surface, can be described by the first 
order multiple regression equation: Z = by + b, X, + b, X,, where bo, b, and b, are partial regres- 
sion coefficients associated with X, and X,. Fitting of curvilinear surfaces is accomplished by poly- 
nomial expansion of the two independent variables in the initial first order equation. 

In the present study, the normalized log-abundances of collembola were fitted by TSA to com- 
binations of spoil moisture content and organic matter. Regression coefficients for first through 
fifth order surfaces were calculated using a FORTRAN program devised by Davis (1973). When 
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several regressions for a particular species yielded significant F-statistics, the equations were tested 
in pairwise fashion in ascending order (i.e. 1st vs 2nd, 2nd vs 3rd, etc.). The residuals about the 
response surface selected as the best fit were then ordinated. 

R-mode principal components analysis (PCA) was implemented prior to, and after TSA, using 
the SPSS subroutine FACTOR (NIE et al. 1975). VARIMAX rotation was applied to the species 
ordinations in order to concentrate as much variation as possible onto the fewest number of principal 
component axes. In each monthly data set, core samples were considered as individual quadrats or 
stands. Factor scores for the stands were secondarily computed from the factor loading matrices for 
the species. The factor scores derived from PCA conducted on the data sets before TSA will be sub- 
sequently referred to as the **pre- TSA" stand ordinations. Factor scores obtained from PCA of the 
matrices of residuals remaining after TSA will be subsequented referred to as the detrended or “post- 
TSA” stand ordinations. For each month, factor scores for the first two principal components of the 
"pre-TSA" and *post-TSA" stand ordinations were rank correlated with water potential, water 
content, organic matter, pH, and predicted sediment density. 

Factor scores for the first four principal components of the **post- TSA" stand ordinations were 
also subjected to Multiple Discriminant Analysis (MDA), with group membership based on vege- 
tation classes. Linear discriminant functions were calculated using the SPSS subroutine DISCRIMI- 
NANT (Nm: et al. 1975). Since equal sample size among groups could not be assumed, the Wilk’s 
lambda likelihood ratio generated for the first discriminant function was used as a crude test of 
equality among multivariate means; the discriminant function coefficients are equivalent to coeffi- 
cients used in multivariate analysis of variance (Morrison 1976). The vegetation classes served as 
"treatment" levels in Kruskal-Wallis tests of the stand scores on the individual components. In 
addition, linear dependence of stand scores on vegetation class was examined by rank correlating 
the two variables. 


4. Results 


4.1. Vegetation classes 


With the exception of pH in June, and water potential in August and September, there 
were significant differences among the vegetation classes with respect to mean levels of the 
other measured environmental variables (Table 1). All the environmental variables included 
in Table 1 were significantly correlated with vegetation class. Vegetation class decreased 
with increasing slope, very strongly so in May (rs = —0.84, p < 0.001, n = 45), and to a 
lesser extent in the ensuing months (Chi square test of homogeneity of coefficients, p — 0.05; 
mean rs = —0.52, p < 0.001, n = 180). Light moss cover (“Im”) tended to be associated 
with high mean slopes (and therefore, high microrelief), independent of plot location. Simi- 
larly, light litter cover ("LL") was restricted to areas of high microrelief, which were sub- 
jected to greatest moisture stress, as evidenced by the low mean water potentials for *LLIm". 
Areas of low microrelief had both thicker litter layers and a more extensive moss turf: these 
areas had higher mean water potentials, soil organic matter and pH levels (Table 1). Pre- 
dicted sediment density also differed significantly among vegetation classes for all months 
(Kruskal-Wallis test, p — 0.001); since sediment density could not be estimated for all 
plots, the significantly different mean rank sums from the non-parametric ANOVAs have 
been arbitrarily denoted “high”, “medium”, “low” and “very low" (see Section 4.3). Rank 
sediment density decreased linearly with increasing vegetation cover (Chi-square test of 
homogeneity of coefficients, p > 0.05; mean rs = —0.69, p < 0.001). 


4.2. Components of spoil water potential 


The lowest (i.e. most negative) water potentials coincided with the lowest spoil moisture 
content measured during the 1980 field season in May. DiLLon (1973), in an earlier study of 
plant germination success at Luscar, observed that a period of soil moisture deficiency 
occurred immediately following April to May snowmelt. Results obtained in the present 
study confirmed his prediction; in mid-May, water potentials in the 0 to 5 em stratum of 
Luscar spoil averaged —1.37 + 0.19 MPa on the 8 year old plot, —5.44 + 0.18 MPa on 
the 4 year old plot, and —4.06 + 0.46 MPa on the 2 year old plot (W. F. J. Parsons, 
personal observation). 

Frequent, light rain likely lessened moisture stress on all three plots towards mid-summer, 
as evidenced by increases in both mean water potential and mean moisture content (Table 2). 
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Table 1. Kruskal-Wallis tests of edaphic variables according to vegetation classes assigned to sample 
cores, May to September 1980 


Vegetation Class!) Kruskal-Wallis 


Month LLIm MLIm MLmm HLIm HLmm K statistics Ts) 
May y —4,31»  —3.39ab — 1.42 1.222 — 174.74 — 0.44 
SOM?) 3.758 3.8748 8.19¢ 6.45» 21.60 0.47 
pu 7.250 6.69 6.484 6.73ab 21.19 0.72 
Slope) 5.3e 3.48 4.1» 2.9 45.35 —0.84 
Sed?) High Med Low Low 17.07 —0.67 
June y — — — — -- -— -- 
SOM 3.644 3.488 4.014 5.43» 5.76 9.71 0.39 
pu 6.94 6.61 6.51 6.51 6.78 7.22 NS —0.43 
Slope 5. 8e 5.56 3. 95 3. 9 2.7 12.90 — 0.51 
Sed High Med High Low Low 32.59 — 0.66 
July y —].43» —0.65ab 0.373» — 1.133» .—0.32* 12.08 —0.50 
SOM 3.902 4,48be 3. 98 b 6.374 5.06¢ 10.05 0.50 
pH 7.06¢ 6.508 6.74 7. 00e 6.492 18.77 —0.58 
Slope 5.9¢ H.1be 4.6» 3.08 3. 3a 15.23 —0.55 
Sed High Med Low Low Low 22.14 —0.67 
August y —1.239  —1.08 — —0.83 . —0.10 434 NS —0.31 
SOM 3.118 3.915 — 4.97¢ 5.184 10.29 0.43 
pH 1.08€ 6.538 — 6.89» 6.77 9.50 —0.39 
Slope 6.1e 45^  — 3.48 3.32 20.57 — 0.62 
Sed High Med Low Low 171.24 — 0.69 
September —0.57 0.50 | —0.23 —0.11  —0.44 5.42 NS —0.35 
SOM 3.58 2.78 3.56 5.50¢ 5.914 11.64 0.34 
pH 7. 03e 6.79» 6.684 6.82 6.768» 12.48 —0.31 
Slope 5. 7e 4.7 4.5 2. 7a 3.62 14.52 —0.42 
Sed High Med V. Low Low V. Low 24.32 —0.75 


1. Upper case letters refer to Light (L), Medium (M), and Heavy (H) litter cover, lower case letters, 
to light (1), medium (m), heavy (h) moss cover, corresponding to designations assigned to stand 
clusters in Figs. 2 and 3. Means in each row superscripted by the same letter do not differ signifi- 
cantly at p = 0.05, as determined from Kruskal-Wallis non-parametric ANOVA and post hoc 
contrasts. K-statistics that do not exceed critical Chi-square values are denoted “NS”. 

2. Spearman rank correlations (rs) between the edaphic variables and vegetation classes are signifi- 

cant at p = 0.05. 

SOM — percentage soil organic matter. . 

. Slope class 

. Predicted sediment density. *Med" — medium, “V. Low" — very low. 


Ot emo 


Typical water potential values of —2.0 MPa to —0.1 MPa were obtained during the period of 
initial spoil moisture deficit. However, possibly as a consequence of collecting samples im- 
mediately following an extensive period of rainfall (late-May to mid-June), 33% of the June 
water potential readings fell outside the reliable range of the instruments employed in the 
determinations. Even after elimination of these measurements from the June data set, the 
regression equation water potentials from percent spoil moisture content was not significant 
at p = 0.05, and therefore, has not been included in Table 2. However, there were significant 
differences between months as indicated by Scheffé tests of the dummy variables included 
in the preliminary regression analysis, and as indicated by the magnitude of the standard 
errors of the monthly regressions, and the partial regression coefficients (Table 2). 

Water potential and water content were significantly and inversely related (product- 
moment correlations, p < 0.001, n = 45) in May (r = —0.73), July (r = —0.74), August 
(r = —0.83), and September (r = —0.55). However, simple linear functions of water con- 
tent, together with additive or multiplicative combinations of the other environmental 
variables, were not sufficient to predict moisture availability (spoil water potential) over 
the wide range of conditions present on Luscar spoil. Trend surface analysis of the water 
potential data was also not useful in this regard. Power functions (as depicted in Table 2) 
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Table 2. Moisture release curves based on of log-transformed psychrometric water potential versus 
log-transformed spoil moisture and organic matter content, May to September 1980 


Month Regression Coefficients (+ SE) Water Water 


Content Potential 
b, b, b, R? F SE» (95) (MPa) 
May 1.45 1.05 0.41 0.77 70.90 0.22 6.60 — 3.63 
(0.09) (0.14) (0.14) (1.10) (0.31) 
June — -- — — — — 15.96 ND 
(1.00) 
July 1.94 —1.76 — 0.82 199.39 0.21 14.23 —1.26 
(0.15) (0.12) (1.09) (0.13) 
August 2.80 — 2.42 — 0.85 233.80 0.22 17.79 — 1.07 
(0.20) (0.16) (1.08) (0.13) 
September 2.41 —2.07 — 0.92 497.25 0.20 36.51 —0.51 
(0.15) (0.09) (1.20) (0.12) 


Linear regressions of 45 replicate measurements per month are significant at p = 0.001. SE» is the 
standard error of the regression. 


Note: 

b,  Y-intercept, the air-entry coefficient of the spoil (logio MPa) at 1% moisture and 1% organic 
matter content. 

b, partial regression coefficient, logy) [0v/0s]. 

b, partial regression coefficient, logy, (soil organic matter). 

ND Water potentials measured in June have been excluded from the regression analysis. 


were the most efficient means of detailing relationships between water potential and the 
other physico-chemical variables; most of the variation in logio y was accounted for by 
either logy) (water content) alone, or logio (water content) plus logo (organic matter). Inclu- 
sion of other log,9-transformed variables did not significantly increase the variance explained 
by regression. Regression analysis using the untransformed variables required estimation 
of 4 or more regression coefficients, and explained a smaller proportion of the variance. The 
power function equations summarized in Table 2 were also consistent with empirical models 
of soil water relations devised by soil physicists (ef. HILLEL 1971). 

Progressive weathering and settling of the Luscar spoil likely altered the spoil moisture 
retention characteristics, and subsequently, the moisture release curves summarized in 
Table 2. There was a significant increase (2-way ANOVA F,.4, = 13.77, p < 0.001) in bulk 
density from a mean 1.17 g em- in May to 1.54 g cm~ in September 1980. Consequently, 
mean available pore space decreased from 50 to 34% of total core volume. The decrease 
in pore space may explain the apparent increase in spoil moisture content through time, 
even though volumetric moisture content (Oy) progressively decreased over the field season 
(W. F. J. Parsons, personal observation). 


4.3. Litter bag burial and microarthropod colonization 


Although spoil erosion was evident on the 8 and 4 year old plots during 1980, there were not 
enough litter bags to estimate sediment density on the two plots using least-squares regres- 
sion. Actual sediment densities measured on the 8 and 4 year old plots were very low, usually 
< 100 mg m-?. Predicted sediment densities on the two older plots were estimated by inter- 
polating between grid locations with known sediment densities. Regression analysis was 
only used on data obtained from the 2 vear old plot; there were no significant differences 
(Scheffé tests, p > 0.05) in mean sediment density among months on the 2 year old plot. 
Actual sediment displacement and redeposition amounted to 6.74t ha-! mo (+ 1.34 
x 103 t ha-! mo~?) as determined from the litter bag measurements averaged over 4 months. 

The depth to which litter bags were buried by spoil was reflected in the moisture content 
of the sediments. Actual sediment density was significantly correlated (p — 0.001, n — 60) 
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with sediment moisture content (rs = 0.85), litter mass losses (rs = 0.70), and the organic 
matter content of the sediment (rs = 0.57) during every month of 1980 when litter bags 
were collected. The commonest microarthropods colonizing the buried litter bags, the 
collembolan Hypogastrura denticulata (BAGNALL), were significantly correlated with 
sediment density (rs = 0.68) and sediment moisture content (rs = 0.67). A second hypo- 
gastrurid, Hypogastrura manubrialis group cf. brevispina (HARVEY), was also significantly 
correlated with sediment density in July and September 1980 (mean rs = 0.61, p < 0.001, 
n = 40). 

Increased collembole densities with increased sediment density reflected an amelioration 
of moisture stress on Luscar spoil. The moisture content of the spoil-infiltrated litter bags 
was 30 to 40% (O. D. W.), four to ten times higher than moisture contents of spoil cores 
taken from the surrounding subplots on the 2 year old plot. Typical water potentials measured 
on the spoil cores collected from the 2 year old plot ranged from —3.6 MPa to —10 kPa 
throughout the field season. Water potentials of the sediments (calculated from pressure 
plate determinations) varied between —0.1 MPa and —15 kPa (W. F. J. Parsons, personal 
observation). 

Germination success of grass and legume seed was enhanced by the placement of litter 
bags on the 2 year old plot, especially where litter bags were buried. Larger numbers of viable 
seedlings germinated immediately under litter bags than from seed sprouting next to those 
bags. The spoil in and under the litter bags served as loci for enhanced plant growth and 
microarthropod colonization since spoil moisture and temperature conditions were likely 
ameliorated. Rooting grasses were found in most litter bags on the 2 year old plot 16 months 
(September 1980) after placement in the field. Although the abundance of H. denticulata 
extracted from these bags was not significantly different from abundances recorded for 
those bags without roots (Mann-Whitney U-test, U = 22, p > 0.20, n = 18), the organic 
matter content of the buried bags containing roots was significantly higher (Mann-Whitney 
U-test, U = 72, p < 0.005, n = 18). The population density of H. denticulata was signifi- 
cantly correlated, in turn, with organic matter content of the bags (rs — 0.58, p — 0.02, 
n = 18). Hypogastrura denticulata and H. manubrialis were probably opportunistic, or 
fugitives species as suggested by the spatio-temporal variation in the litter bag colonization 
data. 

Slope angle (a) and slope length (L) explained 50.5% of the variation in the sediment 
density data(F,,; = 29.10, p < 0.001). The multiple regression equation predicting monthly 
downslope sediment loss and redeposition on the 2 year old plot was: 


logy) SED = 3.96 (+ 0.54) xk Tana + 0.64(+ 0.21) s logy) L 


(The Y-intercept term did not significantly differ from zero, and therefore has not been 
included in the equation.) 

Ranked predicted sediment densities, SED, were consistently (Chi-square test of homo- 
geneity, p > 0.05) and inversely correlated from month to month with soil organic matter 
content (mean rs = —0.53, p < 0.001, n = 225), and vegetation class (mean rs = —0.69, 
p < 0.001, n = 225). Organic matter data derived from analysis of the spoil-infiltrated litter 
bags were too few to be useful in the regression analysis. Vegetation class was excluded from 
regression analysis as well since it varied little from the condition of LLIm (“light litter-light 
moss"). 

Predicted sediment density was significantly rank correlated in May, July, August, and 
September with spoil water content (mean rs = —0.45, p < 0.001, n = 180) and water 
potential (mean rs = —0.53, p < 0.001, n = 180); during June, the correlation between 
sediment density and moisture content was positive and not significant (rs = 0.22, p > 0.05). 
Similarly, pH was not consistently correlated with predicted sediment density. The mean 
correlation between the two variables, averaged for May, July, and August, was rs — —0.54 
(p < 0.001, n = 135). The number of species present in the soil cores was inversely propor- 
tional to predicted sediment density in May, July, August and September (Chi-square test 
of homogeneity, p > 0.05; mean rs = —0.45, p — 0.001, n = 180). The correlation of 
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species number and sediment density was not significant in June. However, total Collembola 
abundance and sediment density were significantly and positively correlated in June (rs 
= 0.3561, p < 0.02, n = 45), but negatively correlated in May (rs = —0.36/0.48, p < 0.002, 
n = 45). 

High correlations with vegetation class, moisture content and organic matter suggested 
that SED was related to the erodibility faetor K in the University Soil Loss Equation 
(WiscuMEIER & Smitu 1978). Predicted sediment density was not a measure of erosion 
hazard or annual soil loss per se, since neither standard test conditions were used nor were 
rainfall intensity, duration and frequency known. Knowledge of these factors is critical in 
calculating erosional losses using the USLE equation; however, the equation developed in 
the present study incorporated two basic elements of the USLE, i.e. slope length and slope 
gradient. 

Present annual rates of erosion were conservatively estimated at 8,088 t km-?, much 
lower than the 17,000—40,000 t km-? y-* predicted by Mckenzie & SrupLick (1979) who 
used the Universal Soil Loss Equation to measure erodibility of mine spoil banks in Ohio. 
Their predictions were based on sites with 30? to 35? slopes that had undergone 20 vears 
of erosion, and where the eroding substrate consisted mainly of sand (64°,, by total sample 
mass). 


4.4. Gradient analysis and ordination 


Fourteen species of Collembola were collected from Luscar mine spoil over 2 years of 
study: Arrhopalites principalis Stacu, Bourletiella sy. Hypogastrura denticulata BAGNALL, 
H. manubrialis near near brevispina (HARvEY), Folsomia nivalis (PACKARD) ssp. duodecimse- 
losa HAMMER, F. quadrioculata (TULLBERG), Isotoma ekamani FJELLBERG, I. ef. mani- 
tobae FJELLBERG, I. cf. tigrina Nicouet, I. viridis BourLet, Metisotoma grandiceps REU- 
TER, Onychiurus subtenuis Folsom, Tullbergia sp. and Entomobrya marginata (TULLBERG). 
The “q-scores” (i.e. direct gradient analysis) derived from the abundance data for these 
species were significantly and negatively correlated with soil organic matter in June, July 


Table 3. Summary of Trend Surface Analyses for collembolan species on Luscar mine spoil, May 
to September 1980 


Species May June July August September 
Hypogastrura 1,110 (27, 6.5) 1,980 (12.5, 4) — = 4,385 (16, 2.5) 
denticulata — 5,175 (26, 7.5) — — 4,380 (33. 6) 
— -— — 4,380 (83, 4) 
3 (0.65) 4 (0.70) — — 5 (0.88) 
Hypogastrura — — — — 210 (0, 0) 
manubrialis — -- — — 1 (0.39) 
Folsomia 2,990 (27, 6.5) 625 (4, 7.5) — -— — 
nivalis — 3,455 (12, 7.5) - -— — 
3 (0.79) 4 (0.77) - E — 
Onychiurus 1,000 (21, 5) 790 (11, 7) 1.625 (20, 6.5) 670 (27, 8) 2,590 (58, 5.5) 
sublenius — 1.795 (23, 7) 4,990 (55.5, 10,955 (46, 10) 1,545 (69, 4.5) 
6.5) 
5 (0.79) 5 (0.83) 5 (0.95) 2 (0.64) 5 (0.90) 
Isotoma ` 200 (17.5, -- 1,140 (28, 1) — — 
cf. manitobae 1.5) 
200 (32.5, 10) — — — — 
5 (0.88) - ? (0.63) — 


For the sake of clarity, the regression coefficients defining the response surfaces are not shown, only 
peak densities (individuals per m?), co-ordinates of the peaks (S, SOM), orders of significant surfaces, 
and the correlation coefficients of the equations are included. Correlations significant at p — 0.05 
are singly underscored, and doubly underscored at p = 0.01. 
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Table 4. Spearman rank correlations of 5 environmental variables with the first and second principal 
component scores of the *pre-TSA" (NP) and "post- TSA" (P) stand ordinations, May to September 
1980 


Month Moisture Organic pH Predicted Water 
content matter sediment potential 
density 
May NP I 0.39 0.38 — 0.42 —0.59 —0.55 
II 0.10 —0.05 —0.28 — 0.09 — 0.03 
P I 0.04 —0.08 —0.45 —10.04 0.06 
lI 0.07 — 0.07 0.05 0.12 — 0.04 
June NP I 0.23 — 0.04 0.37 0.54 — 
Il 0.12 0.48 0.15 —0.25 = 
P I 0.07 0.05 0.13 —0.21 — 
I 0.15 — (0.12 0.26 0.50 ide 
July NP I 0.30 0.22 0.51 — 0.62 —0.29 
li —0,22 0.05 0.26 0.42 —0.20 
P | 0.19 0.06 0.40 —0.51 —0.19 
11  —021 — 0.03 0.26 0.37 — 0.22 
August NP I —0.21 —0.19 0.08 0.35 0.31 
II 0.17 0.28 —0.35 —0.40 —0.09 
P I —o —0.20 0.15 0.43 0.23 
II 0.05 0.18 —0.28 —0.28 0.04 
Septem- NP I 0.37 0.24 0.02 —0.45 —0.34 
ber 1 —ti7 —0.10 0.28 0.43 0.19 
P I —0.03 0.03 0.08 0.31 0.01 
II 0.11 —0.04 0.06 —0.13 —0.07 


Note: Underscored correlations are significant at p = 0.05 


Table 5. Separation of stand scores by vegetation class for the first four “post TSA“ principal com- 
ponents, both individually and as a group 


Month post-TSA'' Principal Components Wilks’ Chi-square 
I T III IV lambda 
May!) K 5.686 7.456 5.213 2.485 0.406 34.30 
NS NS NS NS p < 0.01 
Ts 0.14 0.04 0.09 0.12 
June?) K 19.225 30.456 5.468 5.334 0.355 3.24 
NS NS p < 0.01 
Ts 0.16 —0.67 0.35 0.19 
July?) K 42.041 7.685 4.612 6.058 0.406 36.35 
NS NS NS p < 0.01 
Ts 0.48 —0.32 0.08 0.07 
August?) K 6.566 1.869 3.794 2.464 0.485 28.98 
NS NS NS p < 0.01 
Ts —0.35 0.34 0.16 0.21 
September?) I 2.484 12.957 8.795 2.646 0.319 46.25 
NS NS p < 0.001 
Ts —0.15 0.29 0.29 0.07 


Note: Univariate Kruskal-Wallis tests and rank correlations (rs) are included for each component, 
together with the Wilks’ lambda (and Chi-square approximation) for the first discriminant function 
on all four components. Where there was no significant difference among vegetation classes for scores 
along a particular component, the K tests are denoted NS. Rank correlations significant at p = 0.05 
are underscored. 


1) Kruskal-Wallis K-statistics for stand scores on each principal component have 3 degrees of free- 
dom. The Chi-square approximation of Wilks’ lambda has 12 degrees of freedom. 

2) Kruskal-Wallis K-statistics for stand scores on each principal component have 4 degrees of free- 
dom. The Chi-square approximation of Wilks’ lambda has 16 degrees of freedom. 


Fig. 1. Stand scores derived from Principal Component Analysis of collembolan abundance data 
prior to the use of Trend Surface Analysis (*pre-TSA"), May 1980. Closed circles, triangles, and 
squares represent 8, 4 and 2 year old plot samples respectively. 


and September 1980 (Chi-square test of homogeneity, p > 0.05; mean rs = —0.390, p < 
0.001, n — 135). The remaining environmental variables were not significantly correlated 
(p > 0.05) with the sample scores. The lack of correspondence between the **q-scores" and 
spoil moisture content (and water potential) was likely due either to an incorrect choice of 
initial species weightings or to changes in faunal environment relationships between field 
seasons. The results of the correlation analysis suggested that direct gradient analysis was 
not the most suitable method of analyzing the Collembola data sets, 

The results of the ensuing indirect gradient analysis conducted on the sereened collembolan 
data sets have been summarized in Tables 3, 4, and 5, Figures 1, 2, and 3, and Appendix I. 
Only 5 of 14 species remained after the data screening procedure: Hypogastrura denticulata, 
H. manubrialis, Folsomia nivalis, Onychiurus sublenius, and Isotoma cf. manitobae. 

Generally, changes in densities of the five species with respect to spoil moisture and organic 
matter were not monotonic; however, a significant, planar response surface was generated 
lor H. manubrialis density data collected in September 1980. Differing orders of response 
surfaces were required to detrend the remaining species for the effects of variation in spoil 
moisture and organic matter content (Table 3). Only O. subtenuis displayed significant 
dependence upon the two environmental variables throughout the 1980 field season. The 
response surfaces generated by TSA were characterized by one or more density peaks; the 
predicted densities, together with moisture and organic matter data as coordinates of these 
peaks, have also been included in Table 3. When translated into water potentials using the 
equations provided in Table 2, the moisture-organic matter combinations were equivalent 
to moisture stresses from < —250 MPa to —0.03 MPa. (June spoil moisture contents were 
converted to water potentials using the regression equation derived from the July environ- 
mental data set.) The euedaphic species, O. subtenuis, tended to have higher peak densities 
at higher spoil water potentials than hemiedaphie species such H. denticulata. The predicted 
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Fig. 2. Stand scores derived from Principal Component Analysis of collembolan abundance data 
after the application of Trend Surface Analysis (*post-T5A"), May 1980. Clusters of samples are 
based on original ranks assigned to litter cover and moss cover. 


peak densities, pooled for all species, were significantly correlated (rs = —0.95, p < 0.001, 
n — 22) with predicted water potential values. 

As realistic trends can be extracted from random data, Davis’ (1973) rule-of-thumb 
regarding the magnitude of the product-moment correlation coefficient associated with a 
trend surface equation was adopted as a second criterion (in addition to the F-statistic) for 
selecting significant regressions (cf. Section 3.2). The correlations have been included in 
Table 3. Equations with significant correlation coefficients — 0.30 were rejected as spurious 
relationships, and therefore, do not appear in the table. Coefficients 0.30 to 0.70 suggested 
cautious interpretation of trends and residuals, while equations with correlation coefficients 
> 0.70 were regarded as conforming closely to the original data points (Davis 1973). 

The species ordinations for the screened data sets have been summarized in Appendix I. 
Approximately 849, to 98% cf total variation was accounted for by the first and second 
components in the *pre- TSA" species ordinations, while only 759% to 89% of the variance 
was similarly accounted for in the complimentary analyses conducted on the normalized 
trend surface residuals (i.e. the “post-TSA” ordinations). About 11% to 25% of the variance 
in the *post-TSA" species ordinations was distributed among principal components III and 
IV. In the *pre-TSA" ordinations, H. denticulata and H. manubrialis had high positive load- 
ings on the same component. Folsomia nivalis and O. subtenuis were similarly paired on a 
second component. Isotoma cf. manitobae usually had positive loadings on still a third 
component. There was no consistent grouping of species loadings derived from the “‘post- 
TSA” ordinations. 
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Fig. 3. “Post-TSA” stand ordinations of collembolan abundance data, June, July, August and 
September 1980. Symbols and clusters are the same as those described in Figs. 1 and 2. 
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VARIMAX rotation of the “post-TSA” species ordinations typically required 17 to 
19 iterations to converge upon the terminal orthogonal solutions. About 35 to 45 iterations 
were necessary to achieve comparable rotation of the “pre-TSA” species ordinations. The 
VARIMAX procedure maximizes eigenvalues, i.e. the column sum of squared loadings 
representing each factor; however, variation within species accounted for by each of the 
principal components was not maximized. Variation in the row sums of squared loadings, i.e. 
communality estimates, suggested that the data partitioning procedure used in the present 
study did not behave in a consistent fashion from month to month (Appendix I). However, 
the advantage gained by detrending the species abundance data using TSA can be more 
fully appreciated in a comparison of Figs. 1 and 2. 

Species loadings were used to produce stand factor scores, plotted in the reduced ‘‘species” 
or “community” space (sensu Gauen 1982) as depicted in Figs. 1 and 2, and Fig. 3. Although 
not shown in any figure or table, ordinations produced from the data sets which had not 
undergone the editing procedure were characterized by “horseshoeing’’, i.e, extreme in- 
turning or involution of the endpoint stands. The “pre-TSA” stand ordinations were less 
affected by involution than by curvature or “arching” of the ordination as exemplified by 
the May data set depicted in Figure 1; the companion “post-TSA” ordination (Fig. 2) 
illustrated substantial reduction of the quadratic distortion imposed by the first principal 
component upon the second. Typically, the projections of factor scores into two-component 
space in subsequent monthly “post-TSA” stand ordinations also exhibited apparent reduc- 
tion in quadratic distortion (Fig. 3). However, July and August ordinations were least af- 
fected by quadratic distortion; the first principal components of *pre-TSA" and “post-TSA” 
stand ordinations were highly, and positively correlated (July, rs = 0.78, p < 0.001, n = 45; 
August, rs = 0.96, p < 0.001, n = 45). 

Rank correlations between the environmental variables and stand scores indicated that 
species abundances had been effectively detrended for spoil moisture and organic matter 
(as well as water potential) by the TSA conducted in all months (Table 4). On the other hand, 
predicted sediment densities were significantly rank correlated (p < 0.05) with principal 
component I (PC I) for all months, and principal component II for July through September 
("pre-TSA") ordinations. Predicted sediment densities were correlated with at least one 
principal component for *post-TSA" stand ordinations from June to September (Table 4). 

Vegetation classes have been superimposed upon the first two principal components of 
the *post-TSA" ordinations depicted in Figs. 2 and 3. The ordinated residuals were separated 
into 4 (or 5) significantly different groups along at least one of the components in the reduced 
community space (Table 5, Wilks’ lambda likelihood ratios, p — 0.05). Components III and 
IV were not depicted in the figures. The significant Wilks’ lambda for the first discriminant 
function derived from the May data was likely due to inequality among the group variance- 
covariance matrices (Box's M, approximate PF;5,,,:, = 2.744, p < 0.001). In subsequent 
months, the assumption of homogeneity among group variance-covariances was not as 
seriously violated (Box's M, p > 0.01). On the other hand, vegetation classes were not sig- 
nificantly different (p > 0.05) on every component according to univariate Kruskal-Wallis 
tests of the ranked stand scores, and in spite of the presence of occasional significant linear 
relationships between the two variables (Table 5, Spearman rs, p < 0.05). Vegetation classes 
and stand scores were significantly correlated in June (PC II, III), July (PC I, II), August 
(PC II), and September (PC II, III). 


5. Discussion 


It has been observed that microarthropods are more frequently found in rhizospheres of 
individual plants in nutrient-poor soils, with the suggestion that microarthropod species 
composition is ultimately dependent upon the physiological characteristics of the individual 
plant species themselves (BLACKITu 1974). However, collembolan and acarine survival may 
simply depend on the ameliorating effects that vegetation cover exerts on soil temperature 
and moisture stress. In Antarctic polar desert. JANETSCHEK (1967) found that the relative 
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humidity of bare soil at 1 to 2 em depth was 50% while adjacent moss cushions remained 
at 70% R. H. Appison (1980), using Polar Ordination, concluded that plant physiognomy 
was the primary determinant of collembolan species composition and abundance under 
individual tundra plants regardless of species. Differences in mean soil moisture content 
among the three categories of plants studied were most pronounced on a raised tundra beach 
ridge, the driest of her three study sites (AppisoN 1975). The three plant categories, Dryas, 
saxifrage, and crustose lichens, had very different growth forms. Moisture levels under 
Dryas were intermediate between saxifrage and crustose lichen; seasonal mean water poten- 
tial 10 em under Dryas was —2.0 MPa, with a seasonal low of —2.7 MPa (Appison 1977). 
These values were similar to water potential estimates made in the present study. 

From a study of microarthropod populations associated with the root systems of individual 
pasture grass and weed species, CURRY & GANLEY (1976) also concluded that the relationship 
between the vegetation and the soil fauna should be considered more of a microclimatological 
phenomenon rather than a species-specific association determined by root physiology and 
rhizosphere microbiology. In a similar fashion, GopDARD (1979 a, b) classified sample cores 
on the basis of presence or absence of lichens and moss prior to performing Principal Coordi- 
nate Analysis on the microarthropod data obtained from maritime Antarctic moss-turf 
communities. Although he effectively increased the number of *species" analyzed in the 
ordination by dividing each taxon (including the collembolan Cryptopygus antarcticus) into 
different size classes, the analysis did not reveal clear cut groupings in faunal species related 
to vegetational composition. However, species abundance was affected by the type of plant 
cover from core to core. The “pre-TSA” stand ordinations produced from the Luscar data 
sets bear certain similarities to Gopparp’s Principal Coordinate diagrams, even though the 
species abundance data were not broken down by size and age class in the present study. 

Multiple Discriminant Analysis proved useful in crudely delimiting the niches of the five 
species ordinated in the present study on the basis of vegetation classes. A formal analysis 
of niche breadth and overlap was deemed inappropriate owing to the unequal proportions 
of the resource categories used (i.e. vegetation classes). These coarse habitat measurements 
were more useful in the present study than detailed mierohabitat analysis since the patchy 
distribution of vegetation and the associated soil fauna was reinforced by extreme fluctuations 
in temperature and moisture, and by restricted patterns of litter fall and redistribution 
which often characterize adverse environments. 

Despite the use of Trend Surface Analysis to separate major sources of variation in the 
species abundance data, soil organie matter (as correlated with predicted sediment density) 
and spoil water potential remained important influences on the distribution and abundance 
of the five species. For example spoil pH was correlated with stand scores in May and July 
post-TSA“ ordinations. Soil pH can affect rates of ionic exchange between the soil solution and 
particles. Soil solutions became concentrated as spoil moisture decreased during the field 
season. Despite low specific conductance reported by VissER et al. (1979), dissolved salts 
likely exerted substantial osmotie stress upon the Collembola colonizing Luscar spoil. The 
density of one species, Isotoma ekmani FJELLBERG, was significantly correlated with pH 
(Parsons & PARKINSON 1986), but the osmotic and matric components of spoil water potential 
could not be differentiated using psychrometric methods. Thus, edaphic variables likely 
assumed as great an importance as litter cover in influencing collembolan abundance among 
microsites characterized by abrupt transitions between favourable and less-favoured habitat 
conditions. i 

Lithic material and aspect are the primary factors controlling surface temperatures on 
mine spoils in the Rocky Mountains (Harrison 1974); Harrison (1974) noted a 17 °C dif- 
ference between south-and north facing slopes with 26° slopes, with the former exhibiting 
ground surface temperatures as high as 52°C. In the present study, surface temperature 
maxima of 40°C to 60°C were recorded beneath moderate litter cover (medium litter- 
medium moss") on the 8 year old plot, while maxima ranged from 30 °C to 50 °C on mineral 
soil on the more northerly facing 4 and 2 year old plots (W. F. J. Parsons, unpublished 
data). These values are comparable to surface temperatures recorded on bare spoil heaps 
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in central Pennsylvania (DEELEY X BorpEN 1973) and England (RIcHARDSON 1958), as well 
as clear-cut areas of subalpine forest adjacent to Luscar (Day 1963). 

In a study of seedbed characteristics and germination success on Luscar spoil, DILLON 
(1973) found that even slight microtopographic variation could greatly influence the sub- 
strate temperature as well as moisture stress regimes of adjacent sampling sites. He also 
determined that spoil moisture content of shaded hollows was twice as high as moisture 
levels on exposed rises held at the same moisture tension, due primarily to redeposition of 
fine-textured materials in the depressions from upslope sites. Alteration of spoil moisture 
retention characteristics by redeposition of sediments was also observed in the present study. 

Soil moisture storage, depletion, and recharge depend on the complex interplay of topo- 
graphy, soil porosity and permeability, the distribution and physiological status of the 
vegetation cover, as well as meterological phenomena. As water potential decreases below 
—1.5 MPa, water movement in the vapour phase assumes increasing importance (RosE 
1968); vapor diffusion in very dry soils, such as mine spoil, subject to marked fluctuations 
in surface temperature may lead to great differences in water potential over very short 
distances as cooling or heating proceeds. 

The clusters of stands which emerged from the ordinations, together with results of the 
Multiple Discriminant Analysis, suggested that grouped stands shared similar assemblages 
of species or “component communities” (sensu Roor 1973), despite wide differences in geo- 
graphical location of the individual sample cores around the Luscar dump site. According 
to Roor (1973), a habitat patch possesses unique properties associated with characteristics 
of the local environment which generate, in turn, a predictable subunit of the greater com- 
munity (viz. the component community). The component community responds not only to 
the physical and chemical characteristics of the habitat patch, and to the level and constancy 
of available resources, but also to the durational stability of the patch itself (SourHwoop 
1977). It follows that the magnitude and direction of ecosystem-level transfers of energy and 
nutrients vary among component communities. Extrapolating from Roor’s definition of 
“component communities", Reice (1974) proposed that community level processes, such 
as decomposition, should vary along environmental gradients, i.e. from patch to patch. He 
determined that leaf litter disappearance within streams varied with the type and size of 
habitat patches along a velocity-sediment gradient; leaf pack size was also a very important 
determinant of litter breakdown in lotic systems. WHITFORD ef al. (1980) advanced a similar 
hypothesis for a terrestrial ecosystem, Chihuahuan desert. They found that microarthropod 
abundance was directly proportional to the amount of surface litter present (Sanvos ef al. 
1978) and as a consequence, in direct proportion to the rate of decomposition occurring in 
exposed creosote bush (Larrea tridentata) litter. In contrast with studies of litter colonization 
and decomposition conducted in desert (e.g. FRANCO et al. 1979; VOSSBRINCK etal. 1979; 
Santos ef al. 1981), and mesic forest systems (e.g. STANTON 1979; SEAsTEDT ef al. 1983), 
Collembola were the major colonists of the buried litter bags extracted in the present study; 
97° of these microarthropods consisted of one species, Hypogastrura denticulata. 

In the present study, stand scores derived from the ordinations were significantly and 
positively correlated with vegetation class and predicted sediment density. In addition 
species number was significantly and negatively correlated with predicted sediment density, 
and collembolan abundance in the buried litter bags was positively correlated with actual 
sediment density. As in the Chihuahuan desert, Luscar mine spoil was subject to high surface 
temperatures, sheet surface runoff, wind scouring, and intermittent burial of litter. In addi- 
tion to litter burial, litter was probably redistributed to a certain extent by wind and surface 
runoff; light litter cover tended to be associated with high microrelief. The patchy distribu- 
tion of litter can also be related to the ease with which vegetation is established on mine spoils 
with various angles of repose. Errer (1971, in Ditton 1973) found that germination success 
decreased exponentially as slope angle increased. Above 26? to 34? (i.e. slope class value 
of 8), mine spoil was too unstable to allow sufficient numbers of grass and legume seedlings 
to take root and thus, prevent further erosion. DILLON (1973) speculated that the marked 
reduction in germinant coverage was due to erosion of the seed bed prior to germination, 


106 Pedobiologia 29 (1986) 2 


lack of moisture due to increased runoff, and seedling mortality due to abrasion by runoff 
and colluvium. When continual mechanical weathering was also taken into consideration, 
the correspondence between light moss cover and high slope class values (i.e. 7 to 8) suggested 
that differences in moss cover in the present study were due to instability of the supporting 
substrate as well as low moisture levels. Thus, it was quite likely that new microsites were 
being continually created while old patches were destroyed by ongoing erosional processes. 

Physical disturbance has been considered responsible for maintaining the mosaic character 
of many communities (Levi & Parne 1974; Sovsa 1984). It has been suggested (Lovcks 
1970) that species diversity increases when disturbance maintains competitively inferior 
species (i.e. fugitive species) which survive by colonizing newly disturbed patches. The 
hypothesis has been advanced (CONNELL 1978) that diversity is maximized by intermediate 
levels of disturbance, and that only fugitive species, such as H. denticulata survive under 
high disturbance regimes. 


6. A postscript on ordination in soil faunal research 


Ideally, ordination algorithms should make provisions both for non linear species response 
and for correlations between environmental gradients. While a direct comparison of species 
abundances (criterion variables) against physico-chemical variables (predictor variables) 
would appear to be the most logical procedure, this method of Canonical Correlation Analysis 
(CCA) is not robust to the departures from normality and linearity generally encountered 
in complex communities (Gavcu 1982). At the other extreme, Gaussian Ordination (GO) 
fits normal functions to species abundance data along environmental gradients (Inu & van 
GROENEWOUD 1975). While more realistic in the context of models of community organiza- 
tion, GO cannot be easily extended to more than one gradient when data display substantial 
non-systematic variation (Gauen ef al. 1974). In addition, the normal curves generated by 
GO poorly approximate more skewed and multimodal data. 

The most efficient and biologically meaningful method of ordinating community data 
sets lies between the extremes of common factoring techniques and complex, non-linear 
manipulations such as GO. A suite of techniques, under the umbrella of Catenation (Gauci 
1982), can be brought into play to maximize the analytical power of the established factoring 
methods: distortions in the orthogonal components or.abstract axes can be removed by: 
Parametric Mapping (Noy-Merr 1974), Polynomial Ordination (PHILLIPS 1978), and De- 
trended Correspondence Analysis (HILL & Gauca 1980). CARLETON (1984) has improved 
upon Jounson’s (1981) technique by combining Detrended Correspondence Analysis with 
Canonical Correlation. He has formalized the procedure for simultaneously partialling-out 
unique variation in species abundance data while conducting the Canonical Correlation 
Analysis upon environmental data and ordinated community data sets as Residual Ordina- 
tion Analysis. In this way, he has obviated the problem of correlation between gradients 
that plagued earlier attempts at indirect gradient analysis. Soil biologists should avail them- 
selves of the technique, since certain problems, e.g. microarthropod aggregation, have not 
been satisfactorily explained. Field data suggesting dependence of soil faunal populations 
on water availability, and laboratory experiments demonstrating food preferences may be 
manifestations of the same influence: preferred foods may grow best at a particular field 
water content (UsuER 1976). Detrending procedures could help elucidate relationships be- 
tween the soil microarthropods and the microflora upon which they feed. If microarthropod 
and microbial data matrices were both detrended for some common source of variation, it 
might be possible to summarize competitive or predator-prey relationships within a multi- 
variate space of few dimensions. 
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Appendix I Summary of Principal Components Analyses of screened monthly Collembola data sets, May to September 1980 
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Component loadings and communalities (h*) are given for each species, together with eigenvalues (A) and variance explained by each component (%). 


Note: Species codes = HYDE Hypogastrura denticulata; FONI Folsomia nivalis; ONSU Onychiurus subtenuis; 


Hypogastrura manubrialis 


ISMA Isotoma ef. manitobae; HYMA 
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Synopsis: Original scientific paper 
Parsons, W. F. J., & D. Parkinson, 1985. The use of data-partitioning in the ordination and gra- 
dient analysis of a collembolan community. Pedobiologia 29, 91—111. 

During May to September 1980, 14 species of Collembola were collected from three reclaimed, 
subalpine coal mine spoils differing in age since reclamation (8, 4, and 2 years). Five of the species 
were subsequently used in a community level study employing gradient analysis and ordination 
(Principal Components Analysis). The ordinations were constructed from the screened species ab- 
undance data prior to, and after curvilinear regression analysis (Trend Surface Analysis) on spoil 
moisture and organic matter content. The Trend Surface Analysis was intended to detrend the coll- 
embolan density data with respect to a gradient of soil moisture stress. The residual variation dis- 
played much less pronounced distortion (“arching”) when ordinated in a reduced data space then 
the species abundance data which had not been subjected to regression analysis prior to ordination. 
Factor scores obtained from the ordination of the residuals were related to the density of litter and 
moss cover by Multiple Discriminant Analysis, and were indicative of microsite stability. The factor 
scores were also significantly correlated with a gradient of erosion potential. 

Data-partioning techniques could not completely separate the highly correlated components 
of the complex-environmental gradient of water potential, but Trend Surface Analysis, when com- 
bined with PCA, was a more useful strategy of delimiting habitat utilization in the Collembola than 
the direct application of ordination to the faunal data sets. 

Key words: Collembola, gradient analysis, Principal Components Analysis, Trend Surface Analysis, 
data partitioning, mine spoil reclamation. 
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